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When studying the propagation of VLF and ELF radio waves, Galejs and Row [1964] and Galejs 
[1964] took the ionosphere to be a planar stratified magnetoplasma medium. A wave equation was 
derived which describes the propagation of vertically polarized waves along the magnetic equator. 
It was shown that the wave equation could be solved exactly in closed form for a particular height 
variation of the elements of the permittivity tensor in the ionosphere. The purpose of the present 
paper is to give a method of treating the wave equation which enables other exact closed-form solu- 
tions to be obtained. 



1. Introduction 

In a planar continuously stratified isotropic medium, 
the general electromagnetic field consists of a super- 
position of two independent parts. These correspond 
to horizontally and vertically polarized waves, in which 
the electric and magnetic vectors respectively are 
parallel to the stratifications. The problem of finding 
the field distribution reduces to solving two second 
order ordinary differential equations. A great deal 
of attention has been given in the literature to obtaining 
exact closed-form solutions to these equations for par- 
ticular profiles of the refractive index [e.g., Bremmer, 
1958; Brekhovskikh, 1960; Budden, 1961; Ginzberg, 
1961; Wait, 1962; Gould and Burman, 1964]. 

When a static magnetic field is imposed on an ionized 
medium, the medium becomes anisotropic. In an 
anisotropic stratified medium the electromagnetic 
field cannot, in general, be separated into two inde- 
pendent parts. In this case, the problem of finding 
the field distribution can be reduced to solving a pair 
of coupled second order ordinary differential equations 
[e.g., Budden, 1961]. These equations can also be 
written in the form of four first-order coupled differ- 
ential equations or as a single fourth-order differential 
equation [Budden, 1961]. 

When studying the propagation of ELF radio waves, 
Galejs and Row [1964] considered propagation in a 
planar continuously stratified anisotropic ionosphere. 
The two special cases were dealt with in which the 
waves propagate along and perpendicularly across the 
magnetic equator. In the case of propagation along 
the magnetic equator, the earth's magnetic field is 
transverse to the direction of propagation, as well as 
being parallel to the stratifications in the ionosphere. 
Galejs and Row [1964] showed that, in this case, the 
general field is described by two uncoupled differential 
equations. Horizontally and vertically polarized 
waves then propagate independently. The present 
note will be concerned with this case of propagation 
along the magnetic equator. 



Following Galejs and Row [1964], the #-axis is taken 
to be along the magnetic equator, running from east 
to west. The y-direction is vertically upwards and the 
z-direction is south to north. With the static magnetic 
field in the z-direction the tensor relative permittivity 
in the ionosphere can be written 
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where expressions for €1, €2, and 63 are given, for 
example, by Ratcliffe [1959] and Wait [1962]. For 
a horizontally stratified ionosphere €1, 62, and €3 depend 
on the y coordinate only. For waves propagating along 
the magnetic equator, the fields do not vary with z. 

In the case of vertically polarized waves, the wave's 
magnetic field has a z component, H z , only. This 
can be written 



H z — Y(y) - exp(ik x x — icot) 



(2) 



where k x is a propagation constant, a) is the angular 
frequency, and t is the time. It is found that Y(y) 
satisfies [Galejs and Row, 1964; Galejs, 1964] 
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where k () = co/c, c being the speed of light in free space. 
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Writing 






€1 



(6) 



it is found [Galejs and Row, 1964; Galejs, 1964] that 
f{y) satisfies 



where 
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(8) 



Galejs and Row [1964] showed that (3) could be 
solved exactly in closed form when €i and €2 have the 
following exponential variations: 



and 



€1 = ae°y 



e 2 = be°y 



(9b) 



where a, b, and c are independent of y. In this case 
€1 and e 2 are proportional to the same exponential 
function of height. The purpose of the present note 
is to give a method of obtaining further exact solutions 
to (3). 

2. Method 

The procedure is to choose a form for K(y) so that 
(7) can be solved exactly in closed form. Some suit- 
able forms for K(y) can be found in the work of Rich- 
ards [1959] and Kamke [1948]. These forms will 
contain arbitrary constants. 

The present method reverses the conventional pro- 
cedure which would be to choose the profiles and then 
find the field distribution. Choosing K(y) in this 
manner is equivalent to specifying the functions giving 
the field distribution. The problem is then to find 
profiles of €i and 62 which give rise to the specified 
solutions. Such profiles of €1 and €2 are related by 
the equation 
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where K(y) has the chosen form. Equation (10) is to 
be regarded as a differential equation which is to be 
solved in order to find appropriate profiles of €1 and 
€2. This equation can be integrated as follows. 
First the functions S(y) and T(y) are defined by 
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In terms of these functions €1 and €2 are given by 
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(7) where it is assumed that T# 0. 

When S and T are substituted into the expressions 
for P and Q, (10) becomes 
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(9a) This can be rearranged to give 
=|£|-(Sr)=tf(y) + A°-f- 



+ id^ (logr)+ i 

Hence, provided k x ¥^ 0, 
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s(y)= ^r/ 7, { X(y) + ^~f + l^ (log:r) 



dy 



(log T) 



dy. (17) 



Thus the original problem of solving a rather compli- 
cated second order differential equation has been 
reduced to one of performing a single integration. 

The function T(y) can be taken to be any function 
of y. With K(y) and T(y) chosen, (17) gives S(y) by 
a single integration. Then ei and €2 can be obtained 
from (13) and (14); these expressions will involve the 
arbitrary constants contained in the forms chosen 
for K(y) and T(y). For some choices of K(y) and 
T(y) it will be possible to express the integral for S(y) 
in closed form. 

Any form for K(y) can be chosen for which (7) can 
be solved in terms of known functions. By choosing 
different forms for K(y) and T(y) the method given 
above can be used to obtain exact solutions for the 
field distribution for many different profiles of 61 and 62. 
In general, the profiles of €1 and €2 obtained will depend 
on the propagation constant k x . For waves incident 
from free space onto the anistropic stratified medium, 
k x = ko sin 6, where 8 is the angle of incidence. Then 
61 and e 2 obtained above will depend, in general, on 
the angle of incidence. However, in the examples 
given in the next section it will be found that the 
dependance on k x of the arbitrary parameters involved 
can be chosen so that €1 and e 2 are independent of k x . 
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3. Some Examples 



and 



In this section some particular examples of the 
method outlined above will be treated. The function 
T(y) can be chosen to be any function of y. Here the 
choices T=a, T=ae ay , and T=ay" will be made 



A+By 



a[l + (A+By) 2 ] 



(27b) 



where a, a, and n are independent of both y and k x . which are independent of k. v . For these profiles, the 
In each of these cases the form of K(y) as a function of field distribution is given by (25) where bo is related 
y and k x will be chosen so that S(y) and hence ei and to a and B by (20). 



e-2 will be independent of k, v . 

3.1. Case T=a 

The case in which T=a, where a is a constant 
independent of y and k x , will now be considered. 
Equation (17) now becomes 



S(y) = 



k.r 



K(y) + k%—^\c 



(18) 



Any form for K(y) can be chosen for which (7) can 
be solved in terms of known functions. The additional 
restriction that K(y) can be written 



K(y) = b + b 1 g(y) 



(19) 



will now be imposed. In this equation 6 and 6i are 
independent of y, and g{y) is taken to be independent of 
kj-. Also, b and b\ are taken to depend on k x ac- 
cording to the equations 



and 



bo = k 2 a~ l — k\— ik r B (20) 

b\=— ik, r C (21) 

where B and C are independent of k x . Then (18) gives 



S(y) = f{B + Cg(y)}dy 

which is independent of k x . Then €\ and €• 
independent of k x . 

First example: In this example 



which gives 



giy) = o 

K(y)=bo. 
It is seen from (7) that 

/(y) = exp(±^/ 2 y)- 



(22) 
will be 

(23) 
(24) 

(25) 



Thus the fields are expressed in terms of simple 
exponential functions. 
Equation (22) now gives 



S(y) = A+By 



(26) 



where A is a constant of integration and can be taken 
to be independent of k, r . Hence 
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a[l+(A + By) 2 ] 



(27a) 



If # = 0, then 



I 



a(l+A 2 ) 



€ 2 = 



a(l+A 2 ) 



corresponding to a homogeneous ionosphere. If fur- 
ther A=0, then e\ = a~ l and €2 = 0, corresponding to 
a homogeneous, isotropic ionosphere. 
Second example: In this example 



hich gives 



giy) = y 



K{y) = bo + b x y. 



(28) 



(29) 



When b\ =0, this example reduces to the previous one. 
On writing 



f=-6r 2/3 (*>o+M 



(7) becomes 






U=o 



which is the Airy or Stokes equation. rims 

f(y) = Ai(0, Bi(Q 

where Ai{Q and Bi{Q are Airy functions. 
Equation (22) now gives 



S(y)=A + By+ ]Cf. 



Hence 



ind 



€i = 



a[l+U*-f-£y+iCy 2 ) 2 ] 



A + By+]Cy 2 
: a[l + U + %+|Cy 2 ) 2 ]' 
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(34b) 
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Third example: In this example 

g(y) = y+-ry 2 



(35) 



K(y)=b + b iy +b 2 y 2 (36) 

where 

b 2 =~ik x D y (37) 

D being independent of k x . Thus 

Cg(y) = Cy + Dy 2 . (38) 

When b> = 0, this example reduces to the previous one. 



or 

K(y) = b + bie a v (46) 

where a is independent of y and k x . On writing 

a=— -exphf (47) 



(7) becomes 
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du 2 udu \ u 2 ) 



v 2 = -(4b )la 2 



After writing 



y = ] 



u b\ 

-46 2 ) 1/4 2bi 



equation (7) becomes 



d 2 f^ ( , 1 a 2 \ , . 



(39) 



(40) 



Equation (48) is Bessel's equation. Thus 
f(y)=Z„(u) 



(48) 



(49) 



(50) 



where Z v represents any Bessel function of order v. 
In this case 



S(y)=A + By+{Cla)e«y. 



(51) 



Hence 



where 



€l = 



, 1 

v + - = 



2 (-46: 



TF'V" it). 
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and 



Equation (40) is Weber's equation. Thus [Magnus 
and Oberhettinger, 1949] 



Ay)=D v (u),D- v - 1 (iu) 



(42) 



where the D's are Weber parabolic cylinder functions. 
Equation (22) now gives 



62 a[l + {A + £y + (C/a)e^} 2 ] 
Fi/JA example: In this example 
g(y) = cos (ay) 



Hence 



€i : 



S(y) = ^ + By + 1 Cy 2 + 1 Df. 
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and 



€ 2 : 



^+By-f|cy 2 + i/)^ 



^+£y + |Cy 2 + |Dy 3 
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(44a) 



or 



£(y) = 60 + 61 cos (ay). 



(52a) 



(52b) 



(53) 



(54) 



In this case (7) can be solved in terms of Mathieu func- 
tions. Equation (22) gives 



Her 



6i : 



A+By + ^Cy 2 + \^Df 



(44b) 



S(y) =A+By + (CI a) sin (ay). (55) 



a[l + {A +By + (Cla) sin (ay)} 2 ] (56a) 



and 



Fourth example: In this example 

g(y) = e«y 



(45) 



e 2 : 



A+By+ (CI a) sin (ay) 
"a[l + {4+fly + (C/a) sin (ay)} 2 ]' 



(56b) 
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3.2. Case T(y)^ae^ 
The case in which 

T(y) = tie"" 



He 



f(y) = ZAu) 



(69) 



(57) 



where a and a are independent of both y and k x will 
now be considered. In this case (17) becomes 

S W = X^f «"» fe)"f e-«« + A|+^} dy. (58) 



where Z v represents any Bessel function of order v. 
In this case (63) gives 



S(y) = (B/a) +(A + Cy)e- a v 



(70) 



where A is a constant of integration, which can be 
taken to be independent of k x . Thus 



K(y) can be chosen to be any function of y for which 
(7) can be solved in terms of known functions. Here 
the form of K(y) as a function of y and k x will be re- 
stricted so as to obtain a function S(y) which is inde- 
pendent of k x . Thus, K(y) will be chosen to have the ano - 
form 



€i = 



1 



ae«»[l + {{Blot) +(4 + Cy)e"^} 2 ] (71a) 



K(y) = bo + bie- a » + b»g(y) 



(59) 



ae*y[l + {(B/a) + (A + Cy)e~^} 2 ] ' (71b > 



where g(y) is independent of k x . Also 6 , 6i, and 6 2 When /4=0 and C = 0, the profiles (71) become 

are taken to depend on k x in such a manner as to 
satisfy the equations 



1 



and 



bo = — k\ — ik s B, 
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b x =^-ik v C, 
a 



ik x D 



(60) 



(61) 



and 



ae"v[l+(Bla) 2 ] 



(B/a) 
a&*v[\ + {Bla) 2 ] 



where B, C, and D are independent of k, r . Then (58) 
becomes 



These have the same form as the profiles (9) consid- 
(62) ered by Galejs and Row [1964]. 

Second example: In this example 



= P -2(*y 



S(y) = e~ a " I {Be a " + C + De a "g(y)}dy 



which is independent of k x . 

First example: In this example 



which gives 
On writing 

(7) becomes 

where 



g(y) = o 



K(y) = b + b l e-<">. 



2b\' 2 / ay 



^ + if + (i-4,/= 

da 1 u da \ u 2 ' 



f-=-(4o„)/aA 



(63) 



(64) 



(65) 



(66) 



(67) 



(68) 



After writing 



and 



giy) = e 



K(y) = b + bie-"v + b 2 e- 2c ">. 



2ib'J 2 



/(y) = expf^)-r(a) 



equation (7) gives 



dW 
du 2 
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-1 



4 2ibl' 2 au 
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a 2 4/ u 2 



W=0. 



(72) 



(73) 



(74) 



(75) 



(76) 



This may be compared with Whittaker's confluent 
hypergeometric equation [e.g., Richards, 1959] 

(pw r-o* ok i— 4fi 

dx 2 4 x \x 2 



W=0 



(77) 
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which has linearly independent solutions W ky „,{ax) and 
and W- k , m {— ax). Hence 



f(y) = exp(^Yw ±k , m (±u) 



(78) 



where 



and 



k = 



b, 



2ibl' 2 a 



b 3 =-ik x E (88) 

where B, C, Z), and E are independent of k x . Hence 

S(y) = y -n I {By* + C + Dy n ~ 2 + Ey n g(y)} dy (89) 

which is independent of k x . 
''^' First example: In this example g{y) = and n = 2. 

Thus 



2 ^0 

m 2 = -• 

a 1 



(80) 



Either sign can be taken for m since Wk, +m(x) 
= Wh,-m{x) [Magnus and Oberhettinger, 1949]. 

In the present problem, (63) gives 



On writing 



and 



S(y) = (B/a) + (A + Cy)e~ a v - (/)/a)e~ 2 ^ 
where A is a constant of integration. Thus 

1 

€l ~ ae^[\ + {(Bja) + (4 + Cy)e- a " - (/)/a)e- 2 ^} 2 ] 



(81) 



(82a) 



£(y) = 6o + ^- 



u = bfy 



f(y) = y ll2 Hu) 



equation (7) gives 

d 2 h 1 



du 2 u du 



+ 1--Jh = 



(90) 



(91) 



(92) 



(93) 



and 



where 



€ 2 : 



(B/a) ± (A + Cy)e~ a y - (D/a)e- 2 ^ 
' ae a y[l + {(B/a) + M + Ctfe-"* - (£>/a)e~ 2 ^} 2 ]. 



(82b) 



^ 2 =|" (61 + 62). 



When 62 = and hence D = 0, these profiles reduce to Equation (93) is Bessel's equation. Hence 
those of the previous example. 



3.3. Case T(y) = ay n 
When T{y) = ay r \ (17) becomes 

sw =^M* w -^ + ^h- (83) 

The quantities a and n are taken to be independent 
of both y and k x . When n = the function T(y) con- 
sidered here reduces to T(y) = a, which case was dealt 
with in section 3.1. 

K{y) is now chosen to be of the form 



f{y) = y^Z v (bTy) 

where v is given by (94). 
Equation (89) gives 

S(y) = | By + (C + D)y 1 + Ay~ 2 
where A is a constant of integration. Hence 

1 



ay* 



K(y) = b + ^ + ^ + b :ig (y) 



(84) 



and 



where g(y) is independent of k x . Also 6 , 61, 62, and 
63 are taken to depend on k x through the equations 



e 2 '- 



-\^ByMC^D)y-^Ay- 2 



-By + (C + D)y-' + Ay- 2 



(94) 



(95) 



(96) 



(97a) 



ay- 



bo = — k x — ik x B, 
b\ = kla~ l — ik x C ', 

62 — — .-: rc(rc — 2) — t/cj-Z), 
4 



l+|-fiy+(C+D)y-'+^y- 2 



(97b) 



(85) 

Second example: In this example g(y) = and re=l. 
(86) Thus 



(87) 



K(y) = b + b ^+ b ^ 

y y 



(98) 
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Equation (7) can now be compared with Whittaker's where 
equation (77). Thus 



f(y) = W ±k ..„ l (±2ib}l>y) 



where 



and 



k = 



6, 



2t6J' 2 



(99) 



(100) 



m* = --b 2 . 



Equation (89) gives 

S{y) = \By+C + (A + D\ogy)y-' 
where A is a constant of integration. Hence 

1 



6i : 



and 



ay 



\^By+C + (A-h 



D log y)y x \ 



(102) 



(103a) 



1 



By+C + {A+D\ogy)y- 1 



€ 2 =- 



ay 



1+] -By+C + (A + D\ogy)y- 1 



• (103b) 



Third example: In this example g(y) = and n = — 2. 
Thus 



After writing 



K(y) = b< ) + b,f + 



u=ib\l 2 y 2 






and 



equation (7) gives 



f(y) = y-m h{u) 



(104) 



(105) 



(106) 



du 2 ' 



■1 b 3 + 4i 2 



4 4i6p« 16« 2 



A = 0. (107) 



This may be compared with Whittaker's equation 
(77). Thus 



f(y) = y-«W ±lc , m (±ib\l'Y) 



(108) 



k = K(Mb\' 2 )- 



and 



-=!(!->> 



Equation (89) gives 



(101) Hence 



5(y) = C3^ + ^y 2 -By--Dy- 1 . 



1 



€i=- 



ay 



and 



CyZ + AyZ-By-^Dy- 1 ^ 



Cf + Ay'-By-^Dy- 1 



(109) 



(110) 



1) 



(112a) 



ay 



1 + ( Cy'< + /*y 2 -By-jDy l 



4. Conclusion 



(112b) 



When studying the propagation of ELF and VLF 
radio waves Galejs and Row [1964] and Galejs 
[1964] took the ionosphere to be a continuously 
stratified magnetoplasma with the imposed magnetic 
field parallel to the stratifications. A wave equation 
was derived which describes the propagation of ver- 
tically polarized waves for the case of propagation 
transverse to the imposed magnetic field. This case 
corresponds to waves traveling along the magnetic 
equator. The present paper is concerned with ob- 
taining solutions to the wave equation. 

By the method given in section 2 of this paper, the 
problem of solving the wave equation has been reduced 
to one of performing a single integration. In section 
3 a number of examples have been given showing how 
the method can be used to obtain exact closed-form 
wave functions for various profiles of the elements 
€i and €2 of the permittivity tensor in the magneto- 
plasma. 

Most of the wave functions given in these examples 
were obtained by using the table given by Richards 
[1959, pp. 349-350]. This table gives solutions to 
differential equations of the form y" -\-I(x)y = 0. Fur- 
ther wave functions can easily be found by using the 
same general procedure. 

Finally it will be noted that results for the case of 
an isotropic stratified medium can be obtained from 
the examples given by choosing the arbitrary pa- 
rameters concerned so that the function S(y) = 0. 
Then (13) and (14) give € l = l/T(y) and e 2 = 0. 
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